In order to generate values from this distribution one needs three independent Poisson variables, \(X_1 \sim \text{Po}\left(\lambda_1\right)\), \(X_2 \sim \text{Po}\left(\lambda_2\right)\) and \(X_3 \sim \text{Po}\left(\lambda_3\right)\). Then, \(\left(X,Y\right)=\left(X_1+X3,X_2+X_3\right) \sim \text{BP}\left(\lambda_1,\lambda_2,\lambda_3\right)\).